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ABSTRACT 

We present a survey of bright optical dropout sources in two deep, mult i wavelength 
surveys comprising eleven widely-separated fields, aimed at constraining the galaxy 
luminosity function at z « 7 for sources at 5-10L*(z = 6). Our combined survey area 
is 225arcmin^ to a depth of Jab = 24.2 (Str) and 135arcmin^ to J = 25.3 (4ct). We 
find that infrared data longwards of 2/zm is essential for classifying optical dropout 
sources, and in particular for identifying cool Galactic star contaminants. Our limits 
on the number density of high redshift sources are consistent with current estimates 
of the Lyman break galaxy luminosity function at z = 6. 
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1 INTRODUCTION 

Direct observations of galaxies at very high redshift can 
cast light on the details of structure formation and early 
galaxy evolution. In recent years, the most commonly ob- 
served tracer of star formation at high r edshift has been th e 
Lyman Break galaxy (LBG) population l|Steidel et al.lll999h . 
These sources are selected photometrically based on a spec- 
tral feature interpreted as being the Arcst =1215. 67A Lyman 
break. Neutral hydrogen along the line of sight leads to a 
forest of redshifted Lyman-a absorption lines, blue contin- 
uum flux is suppressed, and a galaxy will exhibit red colours 
in a pair of filters bracketing the break. Depending on the 
relative depth of the imaging in different bands the galaxy 
will be seen to fall in magnitude between the red and blue- 
wards bands, or even to drop below the detection limit. As 
a result, Lyman break galaxies are often termed 'drops' or 
'dropouts' and their approximate redshift is determined by 
the band bluewards of th e break. Hence 'U' - and 'G'-drops 
he at z 3 and z « 4 Jsteidel et al.lll999l), 'V- and 'R'- 
drop s at z « 5 (|Lehnert fc Bremer! 20031; iGiav alisco e t sl\ 
2004 ) and 'I'-drops at z « 6 jStanwav. Bunker, fc McMal3 
2OO3I : iBouwens et aLllioO^ . 

The use of LBG samples to study star formation at 
ever-increasing redshifts has recently encountered a prac- 
tical barrier. At 2 ~ 7 the Lyman break is redshifted to 
1 fim, a spectral region with poor sensitivity in both optical 
CCD detectors and near-infrared instruments. Photometry 
in the relatively broad photometric filters used longwards 
of 1 fim effectively smears out the signatures of sharp-sided 
spectral features such as the Lyman break, making them dif- 



ficult to distinguish from more gradual spectral features such 
as a Balmer break or molecular absorption in low redshift 
sources. The small field of view, high background, low sen- 
sitivity and low multiplexing of current near-infrared spec- 
trographs further complicate the follow-up of Lyman break 
galaxy candidates. 

Sensitive, wide-format near-infrared imagers and multi- 
object spectrographs on the largest telescopes offer potential 
for z — 7 Lyman Break galaxy searches, identifying bright 
near-infrared sources that drop in the optical bands (as dis- 
cussed below in section [2}. Such surveys require knowledge 
of the z = 7 luminosity function in order to estimate the 
required survey depth in a given area. High contamination 
rates in optical-dropout samples also complicate such sur- 
veys. While the most promising candidates will always re- 
quire spectroscopic follow-up, it is clearly necessary to de- 
velop strategies for minimising the number of contaminat- 
ing lower redshift sources from photometry alone. An im- 
portant tool for fulfilling these requirements is the Spitzer 
Space Telescope, which images at infrared wavelengths. How- 
ever the 85cm aperture of Spitzer and l".2xl".2 arcsecond 
pixels of its IRAC instrument lead to blending and confu- 
sion issues at faint limits (mAs > 24). Existing deep field 
surveys incorporate sensitive imaging at wavelengths from 
4OOOA to 8/xm (the four bands of the IRAC instrument). As 
a result, they can provide a first constraint on the surface 
density of bright galaxies at very high redshifts, and a test 
of the ability of Spitzer to distinguish genuine z > 7 sources 
from lower redshift contaminants. 

In this paper we present an analysis of two multiwave- 



2 E. R. Stanway et al. 



length data-sets, identifying and classifying optical-dropout 
sources in eleven widely separated fields, averaging across 
cosmic variance. In section [2] we discuss the colour selection 
criteria which isolate high redshift galaxies. In sect ions [3] and 
|4]we apply such criteria to the ESO Remote Galaxy Survey 
and the Great Observatories Origins Deep Survey respec- 
tively and consider the resulting candidate sources. Finally, 
in section [5] we discuss the implications of this analysis for 
both the 2 = 7 Lyman Break Galaxy luminosity function 
and future wide-area near-infrared survey. 

All magnitudes in this p aper (optical and infrared) 
are quoted in the AB system ()Oke fc Gunn|[l983h . We use 
a flat Universe with Q.a = 0.7, = 0.3 and Ho = 

70/i7okms"^Mpc"^ 



2 COLOUR SELECTION 0¥ Z = 7 GALAXIES 

We employ the Lyman Break technique, identifying sources 
that drop dramatically in flux between the reddest optical 
filter (z-band) and the bluest commonly used filter in the 
infrared (J-band). This criterion is sensitive to galaxies at 
z > 7, at which redshift flux longwards of the Lyman break 
becomes shifted beyond the optical. Given that the mean 
transmitted flux fra ction through the Lyman -g forest a.t z ~ 
6 is less than 5% (jSongaila fc Cowi^ (20021) and may we ll 
fall to zero between z « 6.0 and z « 6.5 l|Fan et al.ll2006h . 
genuine high redshift candidates are expected to have no 
detectable flux in any optical band, as figure [1] illustrates. 

However, other well known populations exhibit similar 
colours, with relatively high number densities, rendering the 
search for high redshift galaxies akin to the search for a nee- 
dle in a haystack. The 4000A Balmer break in old galaxies 
at 2 ~ 2 could potentially be mistaken for a Lyman break 
at 2 > 7, although such galaxies are likely to be detected 
in sufflciently deep optical imaging. Cool late-type stars, 
particularly class L and T dwarfs, are also o ptically faint, 
but p otentially bright in the near-infrared (e.g. lPatten et aP 
l2006l 'l. In order to identify aU or most interlopers we require 
four broadband colour selection criteria in addition to opti- 
cal non-detection. These are illustrated in figures [2] and |3l 
Cuts in the J — K ox J — 3.6/^m will effectively exclude the 
majority of cool T class stars which are blue in these colours. 
Since these two colours are highly correlated, with most of 
the colour difference arising from J — K, strictly only one is 
necessary. However, given the challenge of observing in the 
thermal infrared from the ground, it is often possible to sur- 
vey large areas to deeper limits in 3.6/im with Spitzer/YR.KC 
than in K from the ground, allowing the J-3.6/im colour to 
provide stronger limits in the case of faint or undetected 
objects. 

The true utility of the infrared, however lies at longer 
wavelengths. As figure [3] illustrates, a further requirement 
that sources are flat or slightly red in the Spitzer/YR.AC 
wavebands is likely to exclude some low redshift galaxies 
and the majority of L class stars, which separate neatly from 
the predicted colours of young starburst galaxies at z > 
7. Unfortunately, it is impossible to exclude some late L 
and early T class stars while allowing for photometric and 
intrinsic scatter in the high redshift galaxy population. 

A final source of contamination potentially arises from 
mature (>1 Gyr old), highly reddened galaxies at intermedi- 
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Figure 1. The filters and model spectra considered in this work. 
The filter response curves of the ESO/VLT R-band, Mauna Kea 
Observatories J and K, and Spitzer/IRAC channels, span the 
rest-frame ultraviolet and optical spectra of galaxies at high red- 
shift. We use the models of Maraston (2005), shown here placed 
at z = 7. 

ate redshifts. Such sources have been identified in deep sam- 
ples selected at 3.6/xm and identified as op tical dropouts (e.g. 
iRodigh iero et al. I l2007l : IStern" et al.|[200a ) and have extreme 
colours in both J — K and K — 3.6/xm, as opposed to mature 
galaxies at high redshift (which can have dramatic J — 3.6/im 
colours, but have near-zero colour in J — K. While it is im- 
possible to eliminate such sources from a sample without 
also rejecting highly-reddened, mature galaxies (which show 
a Balmer break) at z > 8, any galaxies with such colours 
must be treated with caution. Given the scarcity of very 
high redshift galaxies, a Bayesian analysis of the redshift 
likelihood distribution is likely to favour the low redshift 
identification. 



3 THE ESO REMOTE GALAXY SURVEY 
3.1 Survey Outline 

The ESO Remote Galaxy Survey (ERGS, Douglas et al, in 
prep) is aimed at identifying 2 > 5 Lyman break galaxies 
in the same dee p imaging fields ima ged by the ESO Distant 
Cluster Survey l|White et al.ll2005l . EDisCS,). Each of ten, 
widely separated fields was surveyed in the optical (VRIZ) 
using the F0RS2 instrument on the Very Large Telescope 
(VLT), in the near-infrared {JKs) using SOFI on the NTT, 
and in the infrared (3.6, 4.5, 5.8 and 8.0/im) using IRAC 
on Spitzer. Typical depths reached were 27.8 in the V-band, 
27.4 in R, 26.6 in /, 25.9 in Z, 24.6 in the J-band and 
23.8 in K {2a, 2" apertures). Typical seeing was 0.6-0.8" in 
the ground-based data. In the IRAC bands, the 2 cr limiting 
depths were measured in a 4.5" aperture and were 25.1 at 
3.6 fim, 25.2 at 4.5 /im, 23.0 at 5.8 /jm and 22.9 at 8.0 /im. For 
full details of the data see IWhite et al. I (|2005l ) and Douglas 
et al (2008). For the purposes of identifying optical dropouts, 
we required full ten band photometric coverage, yielding a 
total survey area of 225.1 arcmin^ (after removing regions 
with high background noise in the near- infrared). Objects 



A Limit on z = 7 LBGs 3 



1.0 



0.5 



-0.5 




-1.0 -0.5 0.0 0.5 

J-3.6 



Figure 2. The near-infrared colours of optical dropout popula- 
tions. The blue lines show the redshift evolution of young star- 
burst galaxies (solid - flat in fy ; dashed — 50 Myr after start 
of burst, dot-dash — 100 Myr). The red line indicates the red- 
shift evolution of an elliptical galaxy formed at 2 = 5 with ex- 
ponentially decaying star formation (r=0.5Gyr). The models of 
iMarasto are used throughout. The purple diamonds in- 

dicat e the measured colours of cool Galactic stars dPatten et al.l 
l200d) . Shaded regions are excluded by our colour selection crite- 
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Figure 3. The colours of optical dropout sources in the 
Spitzer/IRAC wavebands. For explanation of lines and symbols 
please see figure [2] 



were sele cted in the J- band using the SExtractor software 
package (iBertin fc Arno uts 1996) and their fluxes measured 
at the centroid of the J-band detection in each band us- 
ing the IDL 'Aper' routine. For initial detection we required 
that an area of four contiguous pixels exceeded the back- 
ground by at least 1.5 standard deviations, resulting in a 
basic photometry catalogue highly complete for unresolved 
and low surface brightness sources, but also highly contami- 
nated by spurious objects at the faint end. Photometry was 
measured in 2" apertures in the optical and near-infrared 
and in larger 4.5" apertures in the Spitzer/YRAC bands in 



order to optimise signal-to-noise for our unresolved sources. 
The astrometry in the optical images and near-infrared im- 
ages was found to be precise to better than 0.1" (half a 
pixel) and to within 0.6" in the IRAC bands. We find that 
adding an uncertainty of 0.35" (half the seeing disk) to the 
location of the J-band centroid for our compact, low signal- 
to-noise objects has the effect of introducing an additional 
error to the photometry with a standard deviation of 0.04 
magnitudes at the faint limit of our survey. This is negligible 
when added in quadrature with the Poisson noise and statis- 
tical uncertainty in the background at the same magnitude. 

A refined catalogue was constructed of objects with J- 
band detections (>3(7), but no detection at 2 cr in each of 
the V , R, I and Z bands. The V, R and 7-band imaging 
was combined to create a single optical image (reaching a 
limiting magnitude better than 28.5 in each field) and a 
further non-detection required in this image. Magnitudes 
were not corrected for lensing. However, we note that the 
mean weak lensing effect across these fiel ds is to brighten the 
observed magnitudes by 0.1 magnitudes l|Clowe et al.ll20()6h 
compared to their intrinsic magnitudes. Candidate high red- 
shift sources satisfy Jab < 24.2 (3-4(t depending on field) 
after correction for aperture losses. Each candidate was in- 
spected by eye and obvious data artifacts excluded, leaving 
a sample of 35 candidate sources for further inspection. 

By definition our target sources may only be bright in 
a single band, and this measurement is therefore susceptible 
to spurious sources in the detection band. In order to es- 
timate the probability that our optical-dropout sources are 
not merely spurious detections of peaks in the random back- 
ground noise, we repeated the analysis, inverting the J-band 
imaging in the ERGS fields and determining the number of 
negative noise peaks that would have been detected as ob- 
jects to our limit if they were instead positive signals. 

After discounting obviously spurious detections (cosmic 
ray subtraction residuals, regions of extended noise and de- 
pressions in the background in the wings of bright emission 
objects), we find six genuine noise peaks in the 225.1 arcmin^ 
survey area which would have been mistaken for dropout 
sources if positive. Their magnitude distribution, together 
with that of the optical dropout candidates is shown in fig- 
ure U) Making the assumption that depressions in the back- 
ground of this kind are statistically as likely as positive fluc- 
tuations which would reach the photometric catalogue, we 
estimate that a fraction of 17% of our candidate sample, bi- 
ased towards the faint end, may be spurious detections. We 
note that none of the spurious detections are brighter than 
'J — 23.9' (assuming positive flux) and that none of them 
coincide with source detections in any of the other 9 bands 
surveyed. 

Given the faint limits of our data, and the relatively 
broad PSF of the Spitzer Space Telescope in the IRAC bands, 
we also test the probability that random apertures placed 
on the 3.6pim images that form part of the ERGS survey co- 
incide with infrared-detected sources. Placing 300 apertures 
randomly in each of the ten survey fields, we find that 39% of 
photometric apertures are expected to encounter flux from 
contaminating sources equal to a 3 a peak in the background 
level (corresponding to a detection of J = 24.2). 
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Figure 4. The expected magnitude distribution of spurious 
sources due to fluctuations in background noise (hatched), 
compared with the distribution of optical-dropout candidates 
(shaded). The number of spurious sources is based on detection of 
the total number of negative noise peaks in the J-band imaging 
as a function of magnitude. 



3.2 z = 7 Candidates 

This analysis yielded a total of 35 optical dropout sources, 
distributed across 10 equally-sized fields. Of these the most 
densely populated field contained 7 candidate objects and 
the sparsest field only 1 candidate. 

The infrared colours of each candidate were subjected 
to detailed analysis to determine the most likely identifica- 
tion: z > 7 galaxy, lower redshift galaxy or Galactic star. 
Of the 35 dropout sources, 15 were blended at the resolu- 
tion of Spitzer/YRA.C and cannot be studied further, consis- 
tent with the expected 13.7 confused sources given Poisson 
noise in the small number statistics. Deblending was not 
attempted since, at these faint magnitudes, errors in the 
resultant fiux measurement would be dominated by uncer- 
tainty in the colour and isophotal profile of the confusing 
sources. The remaining 20 optical dropouts had clean in- 
frared photometry. Of these, two sources were eliminated 
because of their blue J — K colours {J — K <-0.8), and a 
further eight were rejected as too blue in the J — 3.6/im 
colour ( J — 3.6^m< —1.0). The most likely identification of 
these sources is that they are T-dwarf stars at moderate he- 
liocentric distances (see section (STTJ • Finally, a further nine 
candidates were eliminated from the sample as their colour 
(or the 2 a limit on the colour) satisfied 3.6 — 4.5/im< —0.4. 
These candidates are considered most likely to be either L- 
class dwarfs or old galaxies at intermediate redshift. 

Thus only one candidate of 20 sources studied has 
colours satisfying the colour criteria in figures [2] and |3] The 
colours of this source are given in table [1] and it is shown 
in figure [5] The remaining source heis a magnitude J — 23.6 
suggesting that it is statistically unlikely to be a chance de- 
tection of a spurious noise peak. While the source satisfies 
our formal selection criteria, its colours are close to the blue- 
wards cut-off in 3.6 — 4.5/im colour. This colour is consistent 
with that of an early T-dwarf and, given the photometric 
error bars on this faint source, marginally consistent with 



identification as a 2; > 7 source in some star formation his- 
tories. Hence, in our 225.1 arcmin^ survey area we find an 
upper limit of 1.43 sources, assuming a similar contamina- 
tion fraction for the infrared- blended fraction, lying at 2 > 7 
and a more likely limit of fewer than 1 high redshift galaxy 
to a depth of Jab = 24.2 (3(j). 



4 THE GOODS SURVEY 
4.1 Survey Outline 

The Great Observato ries Origins Deep Survey (GOODS, 
iGiavalisco et al.ll2004l ) is a publicly- available dataset com- 
prising observations from major ground and space-based ob- 
servatories. We make use of the rl.l photometric catalogue 
release which comprises the results of imaging in the F404W 
(6), F606W (v), F814W (i') and F850LP {z') filters of the 
Advanced Camera for Surveys on the Hubble Space Telescope 
{HST/ ACS). We supplement this with the J- and T^s-band 
data obtained as part of the GOODS programme at the 
VLT l|Vanzella et al.ll2006l ) and with Spttzer/IRAC data ob- 
tained as part of the GOODS Legacy Programme (PI: Mark 
Dickinson). As in the analysis in section|3]we measured pho- 
tometry in 2" apertures in the optical and near-IR and in 
4.5" apertures in the infrared. In these large apertures we 
require non-detection to limiting magnitudes of 27.8 in b and 
V, 27.0 in i' and 26.8 in the z'-band and also non-detection 
in visual inspection of the images (since compact sources 
can be detected to more sensitive limits in the high resolu- 
tion space-based data). The near infrared imaging reaches 
2 a limits of J = 25.8 and Ks = 25.2, while the 2 a limiting 
depths were 27.1 at 3.6 /im, 26.6 at 4.5 /im, 24.8 at 5.8 ^im 
and 8.0 fim. 

We restrict our analysis to the GOODS-South field and 
survey a total of 135.0 arcmin^ with full ten-band photo- 
metric coverage. As before, we train our catalogue in the 
J-band and measure the fiux at the near-infrared location 
in Ks and longwards. Optical counterparts are assigned to 
sources if an HST/ ACS detection is identified within 1" of 
the near-IR centroid. Optical dropout candidates were in- 
spected by eye and spurious detections removed. We limit 
our catalogue to sources with a 4ct detection of J < 25.3. 

As in the case of the ERGS survey fields above, we con- 
sider the fraction of confused sources we might expect in the 
Spitzer/IRAC bands through the investigation of randomly 
placed apertures on the 3.6/im image. While the IRAC data 
in the GOODS fields is deeper than in the ERGS fields, the 
GOODS region was chosen for its low foreground contamina- 
tion, while the ERGS survey exploits lower redshift cluster 
fields. Hence the contamination fraction amongst randomly 
placed apertures in the GOODS-S field was a slightly lower 
fraction of 28% coincidence with 3a fiux from other sources. 



4.2 2 = 7 Candidates 

This analysis yielded a total of 14 optical dropout sources. 
Of these, 4 sources (1 with J < 24.2) are blended at the 
resolution of Spitzer/IRAC and cannot be studied further, 
consistent with the 3.9 sources expected from pure con- 
fusion. A further 2 sources are too blue in their infrared 
colours to be high redshift galaxies (i.e. J-3.6^m< —1.0), 
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RA & Declination (J2000) Jab Opt— J^s ^as — Jab Jab — ^As (J ~ 3.6/jm)AB (3.6 — 4.5/^m) 43 (5.8 — 8.0/im)As 

10 54 41.35 -12 44 07.3 23.61±0.26 (3.8) >4.9 >2.1 <0.15 -0.10±0.31 -0.36±0.27 

03 32 25.22 -27 52 30.7 24.73±0.19 (5.2) >3.3 >1.8 1.81±0.22 3.35±0.20 0.45±0.06 -0.14±0.24 

03 32 16.81 -27 49 07.5 24.79±0.22 (4.5) >3.2 >2.1 1.57±0.25 2.99±0.23 0.05±0.06 -0.36±0.24 

Table 1. Photometry of the remaining sources satisfying the colour selection criteria for high redshift {z > 7) galaxies in the ERGS and 
GOODS survey areas. The Opt— J indicates an optical to near-IR colour, measured from the VRI combined image for the first object, 
and from the F435W band in the latter two candidates. The signal to noise of each source in the J-band is shown in brackets. 




Figure 5. Postage stamp images of the three candidate sources listed in table [T] Boxes are 12 arcseconds to a side. In the case of 
1054-1244 the bands shown are V, R, I, Z, J, K, [3.6], [4.5], [5.8], [8.0]. In the remaining two cases V, R, I, Z is replaced by b, v, i' ,z'. While 
all three sources satisfy the colour selection criteria, their behaviour in the infrared is quite different. 



and are tentatively classified as likely T-dwarf stars. Tiie re- 
maining 8 sources are well-detected at 3.6-8.0/im. However, 
while their 3.6-4.5jUm colours are consistent with either high 
or low redshift identification, in six cases (2 of which have 
J < 24.2) their colours at longer wavelengths are rather blue 
(5.8-8.0/im< —0.4) suggesting that a low redshift interloper 
galaxy interpretation is more likely. 

The photometry of the remaining two candidates is 
given in table [T] We note that both sources are detected in 
the J-band at > 5 sigma, suggesting that they are statisti- 
cally unlikely to be spurious detections based on our analysis 
of the ERGS fields discussed above. Both sources are also 
well detected in the K-hand, making a chance coincidence 
of J-band noise with a Spitzer/IKAC source unlikely. Both 
sources show a red spectrum leading to near equal J — K and 
K — 3.6/im colours. As discussed in section [5J this may sug- 
gest either identification as a high redshift candidate or as 
a reddened mature galaxy at intermediate redshifts. Given 
the scarcity of high redshift galaxies, the latter explanation 
is the more likely. Hence, in our survey area of 135 arcmin'^ 
we find an upper limit of 2.8 sources (accounting for the 
fraction of infrared-blended sources) lying at z > 7 and a 
more likely limit of fewer than 1 high redshift galaxy to a 
depth of Jab = 25.3 {4a). 



5 DISCUSSION 

5.1 The Surface Density of L and T-dwarf Stars 

Given the difficulty in unambiguously classifying extremely 
red sources on the basis of photometry alone, it is instructive 
to determine whether the red population identified in this 
work represents an excess over predicted number counts for 
low redshift objects. Specifically, given that early- to mid-T 
dwarfs are expected to be our most significant contaminant 



on th e basis of their observed infrared colours (jPatten et al.l 
l2006h. it is useful to calculate the surface density of T-dwarf 
staro 

The number of T-dwarf stars with spectroscopic 
confirmation in the literatur e remains small (see 
iLooper. Kirkpatrick. fc Burgasserl |2007| . and references 
therein), although this number is increasing rapidly due 
to follow-up of photometrically-identified ca ndidates in 
wide area surveys such as the SDSS (e.g. IChiu et al.l 
|2006| ). iLooper. Kirkpatrick. fc Burgasserl (|2007l ) report a 
lower limit to the surface density of T-dwarfs (T1-T8) 
of 3 X 10"^ pc"^ based on a catalogue of known sources 
within 10 pc of the sun , although they note that they 
may be considerably incomplete due to the non-uniformity 
of the input surveys. A more complete analysis of the 
surface density of T-dwarfs was performed i n the SDSS 
DR1/2MASS overlap region bv lMetchev et all l|2007l ). using 
Monte Carlo simulation to estimate the effect of selection 
functions in both surveys. They estimate a volume density 
of class T0-T8 dwarfs of (7.0 ± 3.1) x 10"^ pc'^ 

Given the typical absolute mag nitude of a class T3 star 
{Jab = 15.3. iMetchev et all lioOTl ). we would be able to 
detect such a star between a distance of 57 pc (closer than 
which our detectors saturate) and 603 pc (for a limit of J = 
24.2) or 1000 pc (for J = 25.3). Given that these distances 
are well within the scale height of the disk, we make the 
assumption that the volume density of T stars is constant 
in this range. From this, we estimate a total enclosed volume 
for early- to mid-T stars of 2224 pc^ in our shallow survey 
and 3796 pc^ in our deep survey area. This corresponds to 



We do not discuss the density of L dwarfs since an essentially 
unknown fraction, varying with distance, of these will be excluded 
by the optical non-detection criterion. 
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a predicted 16±7 T-dwarf stars in our shallow survey and 
27±12 T stars in our deep survey. 

We compare this with the number of sources in our 
dataset classified as likely T-dwarf interlopers on the basis 
of their infrared colours. To a limit of J = 24.2 in the eleven 
survey fields (ERGS+GOODS) we identified a total of ten 
sources as likely T-dwarf interlopers due to very blue colours 
in the Spitzer/YRAC bands. Assuming an equal ratio of T- 
dwarf candidates to z-drops amongst the infrared-confused 
sources, this yields a tentative volume density of (7.9±0.2) x 
10"^ pc-^ 

We note that this is at the high end of the volume den- 
sity estimates discussed above, but is consistent within the 
error bars, given that this is a near-infrared selected sample. 
Given the huge field to field variation in the number density 
of red sources in our shallow survey, and the still signifi- 
cant uncertainty in the surface density of local T dwarfs, 
our results are consistent with at least half of our identified 
z-drops being T-dwarf stars. 

Our results from the deep survey area are less consistent 
with published volume densities for interlopers, with just 
two optical-dropout sources identified as good T-dwarf can- 
didates to a depth of J = 25.3. We note that our deep sur- 
vey is confined to a single region, the GOODS-S field, which 
also appears to be deficient in 2;-band dropout sources at 
J < 24.2 compared to the mean of our ERGS survey fields, 
and that there is a factor of 6 in surface density between our 
richest and sparsest fields. It is possible that the GOODS-S 
sightline probes a void in the local T-dwarf distribution, or 
that the local density is a poor approximation at increas- 
ing heliocentric (and, in this case, galactocentric) distance. 
Given the very small solid angle subtended by our fields, it 
is likely that multiple sightlines through the Galactic disk 
and halo will be required to determine the mean mean vol- 
ume density of such sources as a function of Galactocentric 
coordinate. 

5.2 Implications for the z = 7 Luminosity Function 

While the detection of zero or very few dropout candidates 
to these relatively shallow limits is not entirely surprising, 
it is interesting to consider possible constraints on the lumi- 
nosity function that may be derived from this result. 

The typical luminosity of the Lyman break galaxy pop- 
ulation, M*, is known to evolve with cosmic time, bec oming 
fainter with increasing redshift (iBouwens et al.l 120071 ). The 
highest redshift sample for which the luminosity function has 
been well const r ained is the i'-dropout population at z = 6. 
iBouwens et al.l |2003) found that such a population is de- 
scribed by a Schecter function with M*uv=-20.2±0.2, a = 
— 1.7 and a surface density cf>* = (1-4 ± 0.5) x 10"^ Mpc~^. 

Given no evolution in the shape and comoving volume 
density normalisation luminosity function between z = & 
and z = 7 (914 Myr and 748 Myr after the Big Bang), we 
can predict the expected number of sources at the bright end 
of the luminosity function. Our colour selection (described 
in section \^ is, in theory sensitive over a broad redshift 
range, from z — 6.5 to 2 = 10. However, in an apparent 
magnitude-limited sample, the redshift sensitivity function 
will be biased towards sources at lower redshift, particu- 
larly for a luminosity function with a steep faint-end. Our 
limiting magnitude is 1.6 mag fainter at z = 10 than at 



z = 6.5, a dramatic difi'erence at >L*. We therefore calcu- 
late our survey vo l ume u sing the effective volume method 
after ISteidel et al] l|l999D . accounting for the different lim- 
iting luminosity in each Az = 0.01 redshift shell, using the 
prescription 

Kff(m) = J dzp{m,z)^ 

where p(m, z) is the probability of detecting a galaxy at 
redshift z and apparent 2-band magnitude m, given the lu- 
minosity function and galaxy colour (randomly perturbed 
by an error of 0.3 magnitudes at the faint limit of the sur- 
vey) , and dz ^ is the comoving volume per unit solid angle 
in a slice dz at redshift z. 

Given the p arameters of the z = 5.9 luminosity function 
discussed above l|Bouwens et al.ir2007l ). we find that our shal- 
low survey is sensitive to 8.0% of the total volume enclosed 
between z = 6.5 and z = 10.0 per unit area to the limiting 
luminosity of the faintest object in the sample, while the 
deep survey selection is sensitive to 16.5% of the same vol- 
ume. We note that if we were to adopt instead t he tentative 
value of M* proposed by IBouwens et al.l |20o3) for z — 7, 
we would survey an effective volume only 4.5% and 9% of 
the total volume enclosed by the redshift range sensitivity 
range. 

In our shallow sample reaching J — 24.2, we survey 
an effective volume of 1.0 x 10^ Mpc'^ and reach a typical 
limiting magnitude (at z = 7) of AIuv~-22.7. Assuming 
no evolution in the high redshift LBG luminosity function, 
we are sensitive to galaxies brighter than 10 M* (given the 
z—6 value for M*). Our detection of no good candidates in 
the combined GOODS-f ERGS sample is consistent with the 
predicted 0.001 galaxies. 

While this is not unexpected it is interesting to note 
that recent evidence from z > 5 galaxy populations have 
suggested that the the universe at earlier times was a more 
vigorously active place. The evolution in size of Lyman- 
break galaxies, together with their more slowly evolving star 
formation rates implies that star formation density rises with 
increasing redshift ( Verma et _ al., 2007 ) • T he identification o f 
old stellar populations in z = 6 galaxies (|Evles et al.ll2005l ) 
also implies that starburst activity must have been rapid 
and intense at early times. Hypothetical luminosity func- 
tions for primordial starbursts (which have a top-heavy ini- 
tial mass function and a higher light to mass ratio than more 
conventional stellar populations) imply that a factor of 100 
or more increase might be seen in th e number density o f 
ultraviolet- luminous sources at z > 7 l|Barton et al.|[20M ). 
The number density of z = 7 candidates in our combined 
ERGS-I-GOODS survey area provides a firm upper limit on 
the bright end of the luminosity function, confirming that it 
inhabits the expected region of parameter space. 

Our deeper sample reaches J = 25.3 or a typical lim- 
iting magnitude (at z = 7) of Mj7v'=-21.6. Assuming no 
evolution in the high redshift LBG luminosity function, we 
are sensitive to galaxies brighter than 3.6 M*. In our effec- 
tive survey volume, 7.7 x 10^ Mpc^, we might expect to ob- 
serve 2.4 galaxies, consistent with our detection of up to two 
good candidates. This is sufficient to rule out the most opti- 
mist ic predictions for m etal-free stellar populations at z > 7 
fsee lBarton et al.ll2004l ). The upper limits derived from our 
analysis are consistent with the measured LBG luminosity 
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function aX z = and suggest that the volume density of lu- 
minous (brighter than M*) galaxies remains approximately 
constant (if both our candidates are real z > 7 galaxies) or 
decreases (if one or more is not real) with increasing redshift 
rat her than increas i ng. 

iBouwens et alj l|2004l ) estimated the surface density of 
optical dropout sources in the Hubble Ultra Deep Field and 
identified four likely candidates to a depth of Hab ~ 27.0 
(as opposed to 8 predicted z = 7 galaxies based on the i'- 
drop luminosity function). Given a spectrum fiat in (as 
appropriate for young starbursts), this is equivalent to two 
magnitudes deeper than our limit in the GOODS fields and 
the area of their survey was only 4% of that discussed here. 
This would imply an expected 0.5 z — 7 galaxies meeting our 
deep survey criteria. While we determine an upper limit of 
2 galaxies at « > 7 in our deep survey, our more realistic ex- 
pectation is that neither is likely to lie at high redshift. The 
confirmation of one or zero sources would confirm that the 
HUDF is not significantly under-abundant compared with 
other sightlines at the same redshift. 

Given the varied colour selections employed by differ- 
ent authors, fair comparisons are difficult. We note that 
we would have rejected the majority of massive z > 
4 galaxy candidates proposed in t he GOODS-South by 
iDunlop. Cirasuolo. fc McLurd l|2007l ) due to the presence of 
fiux in the observed-frame optical or as below the reliable 
de tection limit in the J-ba nd. By contrast the requirement 
of lRodighiero et al.1 (|2007| ) that their candidates be fainter 
than Ks = 23.5 would have excluded a number of the candi- 
date galaxies presented here. Their selection method, based 
on a fiux-limited 3.6/im sample, may effectively select mas- 
sive galaxies at these redshifts, and yet is likely to suffer 
incompleteness due to both confusion and insensitivity to 
young, starbursts (which are highly stochastic and domi- 
nate the number counts in h igh redshift Lyman break galaxy 
samples IVerma et al.ll2007l ). Hence the optical-dropout sur- 
vey presented here is complementary to other z > 7 surveys. 

5.3 Implications for Future Surveys 

In the next few years, the UKIDSS Ultra Deep Survey (UDS) 
will survey an area of 0.8 deg^ to Kab ~ 25 and will incor- 
porate a search for high redshift galaxies. Our deep survey 
probes 12.5% of this area and find at most 2 candidates 
which might require spectroscopic follow-up, consistant with 
16 potential z — 7 candidates in the full volume of the UDS 
- a feasible number given a reasonable investment of large 
telescope time. This modest number of galaxies is largely in- 
dependent of evolution in the faint-end slope of the luminos- 
ity function, since the selection is truncated well above M*. 
Given the z = 6 luminosity function, perhaps seven of the 
candidate galaxies identified in the UDS might be expected 
to lie at high redshift, although continued evolution in M* 
with increasing redshift potentially reduces this number still 
further. It is clear that, even with exceptionally deep aux- 
iliary data bluewards and redwards of the selection bands, 
the ratio of candidates to genuine galaxies in a photometric 
sample is likely to be at least 2.3:1. 

However, as the analysis in this paper makes clear, the 
number of contaminant sources in samples derived from 
optical/near- infrared photometry alone is likely to be much 
higher. Depth-matched imaging in the infrared wavebands 



of Spitzer/IRAC are essential for reducing the number of 
near-infrared photometric candidates (which in the UKIDSS 
UDS could exceed a hundred) to a more practical number 
for expensive follow-up observations. 

Wide-field infrared imagers currently being commis- 
sioned, such as Hawk-I on the VLT, are likely to revolu- 
tionise this field, making z > 7 as accessible to observers 
as z < 6 at the current time. However, the fundamental re- 
quirement that auxiliary observations at longer wavelengths 
are obtained above the atmosphere will limit their utility for 
high-redshift surveys after the inevitable decline of IRAC on 
the Spitzer Space Telescope and until its eventual replace- 
ment by the James Webb Space Telescope. 

One prospect for the hiatus period between these two 
events is to extend the lifeti me of Spitzer in a 'warm' 
mode l|van Dokkum et al.ll2007h . While sensitivity would be 
lost longwards of 5 microns, the two shortest Spitzer/IRAC 
wavebands provide sufficient spectral leverage to identify the 
majority of potential contaminants (>80% in this study), 
and in particular to isolate cool Galactic stars from the 
galaxy locus. Inevitably the size of the 'haystack' will in- 
crease with the loss of the 5.8 and 8.0/im wavebands. How- 
ever, the use of a warm Spitzer for surveys would provide an 
archive that would be of use for many years, complementing 
both existing and future, deep ground-based surveys. 



6 CONCLUSIONS 

The key points presented in this paper can be summarised 
as follows: 

i) We have identified an analysed a sample of z-band 
dropout sources in deep, multi-wavelength imaging, consid- 
ering their suitability as potential z > 7 galaxy candidates. 

ii) In an area of 360 arcmin^ (the ERGS-I-GOODS sur- 
veys) , we determine a strict upper limit of 1.4 potential z = 7 
galaxies and a likely limit of zero candidates to J = 24.2. 

iii) In an area of 135arcmin'^ (the GOODS survey), we 
determine a strict upper limit of 2.8 potential z = 7 galaxies 
and a likely limit of no good candidates to J = 25.3. 

iv) These results are consistent with current estimates 
of the luminosity function of Lyman break galaxies at z = 6, 
for galaxies with L > 5L*. 

v) The use of deep data longward of the Tf-band is 
extremely valuable if the number of high redshift candidates 
from deep surveys such as UKIDSS UDS from unmanageable 
to reasonable numbers for spectroscopic follow-up. 
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